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Ian A. Chapman 
Abstract 

We compare two competing theories regarding finite temperature 
wave-function corrections for the process H e'^e~ and for n + v ^ 
p + e~ and related processes of interest for primordial nucleosynthesis. 
Although the two methods are distinct (as shown in H ^ e"*“e“) they 
yield the same finite temperature correction for all n —> p and p ^ n 
processes. Both methods yield an increase in the He/H ratio of .01% 
due to finite temperature renormalization rather than a decrease of 
.16% as previously predicted. 


1 Introduction 

To calculate finite temperature corrections to physical processes beyond the 
tree level, one needs to have a consistent and correct method of temper¬ 
ature dependent renormalization. Unfortunately, there seems to be more 
than one way to do this. Although everyone seems to agree on how to handle 
temperature dependent phase space, temperature dependent vertex correc¬ 
tions, effective mass corrections, and temperature dependent bremsstrahlung 
corrections, they do not agree on how one handles wave function renormal¬ 
ization. 

In this paper, we will examine two competing models. Each model’s effect 
on the ’effective’ projection operator (which in turn will tell us how to han¬ 
dle wave function renormalization) is examined for the general case. These 
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general results are then applied to if —> e’''e“ since this is a fairly simple pro¬ 
cess to understand. We also apply both methods to processes of interest for 
primordial nucleosynthesis, specihcally n + v p + e~ and related processes. 
These processes were hrst analyzed by Dicus,et al.[|l|, and independently by 
Cambier, Primack, and Sherp|. Since these papers were published, however, 
it has been claimed that neither paper handled the hnite temperature wave 
function renormalization scheme correctly. Both Donoghue and Holstein |^], 
and Sawyer 1^ propose models that claim to correct this difficulty. Both agree 
that one must use finite temperature Dirac spinors (and thus projection op¬ 
erators), but otherwise their models seem to differ significantly. 

Before we begin a few comments about gauges are in order. The correction 
of interest will be the wavefunction renormalization part of the width. Since 
this is a gauge dependent quantity, we might wonder if our choice of gauge 
will make a difference in our results. We solve this problem in two ways. For 
most of this paper, we will avoid a specific choice of gauge. Instead we will 
discuss the differences between the models in the most general way possible. 
When we calculate a specific result for nucleosynthesis, this is no longer a 
viable option. In Appendix I, however, we prove that we can talk about the 
differences in the wavefunction renormalization parts in the width in a gauge 
invariant manner. We take advantage of this in Eq.(14). 


2 Generalized Analysis 


As discussed in the previous literature, the presence of a background thermal 
bath breaks Lorentz invariance. At T=0 we can identify the wavefunction 
renormalization constant Z 2 for fermions with momentum p as I-Z 2 = ^ 

Op j/=m 

where S is the self energy. Since Lorentz invariance is broken when Tt^O, it is 
claimed that this identification of a momentum independent renormalization 
constant is no longer correct. Donoghue and Holstein advise us to look at 
the modified propagator instead and to modify the Dirac spinors accordingly. 
Thus the modified propagator becomes 

^ . ( 1 ) 


G = 


m 


where the fermion self-energy, E is 

E = Apo7° - Bpa^ - C. 


( 2 ) 
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A,B, and C are momentum dependent scalar functions. We identify an 
’’effective” momentum, where = ((1 — A{p))p^, (1 — B{p))p^) and an 
effective mass fh where fh = m — C{p). We are instructed to replace the 
fermion momenta and mass with the ’’effective” fermion momenta and mass 
everywhere inside the Dirac spinors. When this is done, and the projec¬ 
tion operator is expanded about the physical mass shell energy, we get the 
following ’’effective” projection operator 


A 


+ _ 
eff - 


p + m 

2po 


( 3 ) 


Converting this into the mass shell p^, expanding about the mass-shell 
energy, and writing this in terms of A,B, and C, we get (noting that Up = 
or the mass shell energy) 


/ 

9((l-2A)p2_(i_2S)p2-(l-^)m2) 

\ 

1 - 

V 

dpo 


2up(l - A) 

pQ=UJpj 


Simplifying this, we hnd 


/ p + m\ 

V 2a;p ) 


( 4 ) 


A 


+ _ 
eff - 



dA 

dpo 


UJp 


P0=uJp 


dB 

dpo 




PO — 


UJri 


dC 

dpo 


m 


PQ=ijJp Pj 


^ + m 

2ujp 


( 5 ) 


Thus the effective projection operator is equal to the original T=0 projec¬ 
tion operator multiplied by a single scalar momentum dependent function. 
Donoghue and Holstein identify this function as 1-”Z2” where ’’Z 2 ” is the 
effective temperature dependent wave-renormalization factor. 

Sawyer gives us a differing way to handle hnite temperature wave-function 
renormalization. Sawyer tells us to start with the same modihed propagator 
as shown in Eq.(l), but then reminds us that the propagator can also be 
written in terms of the effective Dirac spinors as shown by Weldon]^. Thus 

G = — 1 _ finite terms (6) 

Po 

where w is the temperature dependent renormalized Dirac spinor. Using 
this fact. Sawyer is able to perturbatively derive the temperature dependent 
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renormalized Dirac spinors, tc, in terms of the unrenormalized zero temper¬ 
ature Dirac spinors, u. He was then able to use these new spinors to directly 
dehne the effective projection operators. If we follow this procedure we get 


A+ff = —( 7 ° — 




PO — 


as 

dpo 


7 


(7) 


PO — 


When we write Sawyer’s effective projection operator in terms of A,B,C, and 
Up as dehned previously, we get (evaluating po at Up). 



. an ^ 

A + —po 
y dpo , 


f i) + m\ 

V 2cc;p ) 


'dB , an 2 n dC dC 

Po - 7 - T^Po - 

opo opo opo opo 



( 8 ) 


It is immediately evident that the above expression for is in general 
different from the expression for in Eq.(5). Specihcally, we note that 
A^jj is not equal to the unrenormalized projection operator times some scalar 
functional factor. We see this explicitly in the process H e^e~. 


3 H e+e comparison 


The tree-level matrix element for this process in momentum space is 


Mx — 


igm 


-u(e ]v[e 


w 


( 9 ) 


where m is the mass of the fermion, Mw is the mass of the u and 
V are the Dirac Spinors at T=0, and Mt is the tree level matrix element 
for this process. The modihed width considering only those terms due to 
the temperature dependent wavefunction renormalization can be found (to 
hrst order in a) by cross-multiplying the tree-level matrix elements using the 
modified projection operators, A^jj and Ajjy, in place of the normal zero 
temperature unrenormalized projection operators. One then integrates over 
the correct temperature dependent phase space. Since this later part is not 
in dispute, it is easier to simply compare the modihed spin-summed squared 
self-energy matrix elements in both models. (As I mentioned before, I will 
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neglect the temperature dependent mass terms here.) We assume for this 
process that the Higgs is at rest with respect to the background thermal 
bath, and that the Higgs is signihcantly more massive than the background 
temperature. When we use the modihed projection operator for Donoghue 
and Holstein’s model, we get 


_ dB _ dC 2m \ 

4M^ V ^ deoeo deoMn] 

When we use the modified projection operator as described by Sawyer, how¬ 
ever, we get 




\M\^ 


4 ^ 


2M| 



^ dA \ dC 

A + —eo -h -^mMu 
, ocq I oeo 


( 11 ) 


These matrix elements squared yield differing widths when integrated over 
phase space. Thus the two models are distinct for H —> 


4 n -\- u ^ p + e comparison 

Finite temperature corrections have been of considerable interest for n ^ p 
and p ^ n reactions because of nucleosynthesis. Specihcally, one would like 
to know how hnite temperature corrections affect abundance ratios such as 
H/He in the early universe. Since all n —p and p ^ n processes are related 
by virtually identical matrix elements, a comparison of n + u ^ p + e~ is 
sufficient for all such processes. 

In this comparison we make the following standard assumptions (valid 
for proton to neutron ratios in the early universe). First we assume that 
the ambient temperatures of interest are of order 1 MeV so that the proton 
and neutron have considerably more mass than the ambient temperature. 
(Thus any terms of the form or ^ are ignored.) We also assume that 
the initial nucleon is at rest with respect to the surrounding medium. The 
neutrino electron part of the spin-summed squared matrix element looks like 

T'' (12) 

In general has three kinds of terms. There are scalar terms, terms 
that are multiplied by a factor of 7 °, and the remaining terms which are all 
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multiplied by a factor of 6*7*. The scalar terms will trace away, and those 
terms that are multiplied by do not survive the phase space integration. 
Thus the only terms in the effective projection operators that matter are 
those terms which are multiplied by 7'^. We now look back at Donoghue 
and Holstein’s effective projection operator (Eq.( 5 )) and Sawyer’s effective 
projection operator (Eq.(8)). Isolating the 7° terms in both models, we get 


AJ/ = 


^ ^ 


dB el 


dC m' 


deo Co deo Cq ^ 


eo7 


(13) 


i.e., the effective projection operators are identical. In the above eqnation, 
is the magnitude of the 3 -momentum of the electron. To verify this, I did 
this calcnlation using Donoghne and Holstein’s method and Sawyer’s method 
using the Feynman gauge. (Sawyer used the Coulomb gauge in his analysis.) 

When we calculate the width according to Sawyer (and Donahne and 
Holstein) minus the width according to Dicus et ah (and Cambier et ah) for 
the reaction n + n ^ p + e~ we get 

^■P Gl{gl, + 3gl)a f ^ f ^ 

OTSawyer “ OTDicus = - 7^ - / [Q - e) e„de / k dk 

Zti Jq Jo 


[1 -I- exp{—e/T)] 


-1 


1 -f- exp{ 


e-Q, 
T ' 


-1 


X 


nf{Ek 

Eu 


4m^(e -T E^) 
{E^e + w?y — elk"^ 


Am?{e — Ek) 


Eh 

j_ Elii 

{Eke — rn^Y — elk'^ e^k 


{Eke — eykY — wd 


{Eke + e^aky - 


e ^ 

1_ T y-| 

{Ele^ — {eyk — nYY 

1 nf,(A;)Ln 

e+e„ 

e—Ev 

1 —ijli 

e^k 

E^e^ — {eyk + m?Y 

J eyk 


(14) 


where Q is the rest mass of the neutron minus the rest mass of the proton, 
- m?, Ek = Yk‘^ + w?, nf{Ek) = [1 + exp{Ek/T)]~^ , nh{k) = 
[exp{k/T) — l]“b Using the above eqnation and its analogs for the other 
hve n ^ p and p ^ n reactions, we obtain the change in the width ratios in 
hgures one and two. 

In hgure three we combine the n ^ p and p ^ n processes into a single 
ratio change for n ^ p and p ^ n. 
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We use these results to solve for the change in the asymptotic neutron to 
baryon ratio, 6Xn as described by Bernstein[^] where 

5X,,{t)=S{t) (15) 

where S{t) is the solution of ^ + Tp^n)S = 0. 

We solve Eq.(15) for the 5Xn change we expect. We solve for S'(f), 
and STn^p using Vegas (a Monte-Carlo numeric integration routine). We 
solve the outer integration over t by Gaussian quadratures. We obtain an 
increase of 6Xn of 0.000415. Thus in Dicus et ah’s original paper, Y 3 — W 
in Table 1 (on page 2701) should be changed from -1-0.0004 to -0.000015. 
This corresponds with an increase of .01% in SX^ due to hnite temperature 
renormalization as compared to a .16% decrease as predicted by Dicus et 
ah and Cambier et ah Sawyer predicts an increase of .02% in 6Xn due 
to hnite temperature renormalization. Although his result differs somewhat 
from ours, we both agree that the change in 6Xn should be very small and 
positive i.e., increase due to temperature dependent renormalization effects. 


5 Summary 

It is clear that the models proposed by Donoghue and Holstein, and Sawyer 
are indeed distinct. However, it is equally clear that both models yield the 
same hnite temperature corrections for n ^ p and p ^ n processes when the 
temperature is much less than the nucleon mass. This is because the later 
processes, at least as described in the early universe, are insensitive to any 
correction that does not depend directly on the 7 ° or energy component of 
the ehective projection operator. It is extremely interesting to note that the 
ehective projection operators for both models are identical in this one crucial 
component. 

If one can show that in one’s process of interest, the 7 * and scalar pieces of 
the renormalized projection operator don’t matter, then you can indeed use 
Donoghue and Holstein’s ehective ’’Z 2 ” in the conventional way. Otherwise, 
it is not clear which method, if either, should be used. 
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7 Appendix I 


We want a gauge independent way of comparing two hnite temperature renor¬ 
malization models. It is not clear that looking solely at the wavefunction or 
self-energy part of the renormalized width will be sufficient since this part 
is gauge dependent. Fortunately we know that all of the overall models we 
consider in this paper are gauge invariant and that the only part of these 
models that differ are the self-energy parts. Thus 

^-ptotal = XrSE , ^-pother /ij-N 

Sawyer Sawyer ' Sawyer \-^^J 

and 

^-ptotal _ ^-pSE rpother M V'l 

Dicus — Dicus “T Dicus AC 

In the above equations, represents the total change in width due to 

hnite temperature renormalization, represents the change in width due 
to hnite temperature wavefunction renormalization, while represents 

all other changes in the width due to hnite temperature renormalization. 
Since the models proposed by Dicus and Sawyer only disagree on the wave- 
function renormalization part, we conclude that must be the same for 

both models. Thus if we know subtract Eq.(17) from Eq.(16) we get 

( 18 ) 


^-ptotal _ ^-ptotal 

Sawyer Dicus 


_ xySE _ xySE 

Sawyer Dicus 


Since we know the total change in widths is gauge independent for both mod¬ 
els, it follows from Eq.(18) that ST^sLyer-^'^mcus must be gauge independent 
as well. This gives us a good way to compare two diherent models for hnite 
temperature renormalization. 
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Figure 1: These are the ratios of the change of the hnite temperature hrst 
order corrections (as calculated in Eq.(14)) to the hnite temperature un¬ 
corrected width for aX\ n p processes in nnits of The temperature is 
measured in nnits of the electron mass. The short dashed line represents 
n peu, the long dashed line represents ne pu, and the solid line repre¬ 
sents nu pe. 10 
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Figure 2; These are the ratios of the change of the hnite temperature first 
order corrections (as calculated in Eq.(14)) to the hnite temperature uncor¬ 
rected widths for all p —n processes in units of The temperature is 
measured in units of the electron mass. The short dashed line represents 
pu ne, and the long dashed line represents pe nv 
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Figure 3: These are the combined changes in the hnite temperature hrst order 
corrections (as calculated in Eq.(14)) to the combined hnite temperature 
uncorrected widths for both n ^ p and p ^ nin units of The temperature 
is measured in units of the electron mass. The solid line represents the 
combined n ^ p processes while the dashed line represents the combined 
p ^ n processes. 12 



